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I. INTRODUCTION 



At first glance, Yang-Mills (YM) theory and gravity seem to have little in common. This 
is especially true from the point of view of calculating perturbative scattering amplitudes: 
then the Einstein-Hilbert Lagrangian contains an infinite number of vertices, but the Yang- 
Mills Lagrangian contains only three and four-point vertices. While it is often said that both 
theories are gauge theories, the enormous differences between the theories from the point of 
view of computation makes this point of comparison seem rather poetic. 

But a deeper study reveals concrete relationships between gauge and gravity scattering 
amplitudes. Early in the first heyday of string theory, Kawai, Lewellen and Tye (KLT) 
derived a relation, valid in any number of dimensions, expressing any closed string tree 
amplitude in terms of a sum of products of two open string tree amplitudes [Ij. Taking 
the field theory limit of the string amplitudes, these KLT relations imply in particular 
that any graviton scattering amplitude can be expressed in terms of a sum of products of 
colour-ordered gluon scattering amplitudes. The KLT relations are very simple for low point 
amplitudes; for example, the three-point graviton scattering amplitude is simply the square 
of the three-point Yang-Mills amplitude.^ For higher-point amplitudes, however, the sum 
in the KLT relations becomes complicated, presenting an algorithmic obstacle to computing 
graviton amplitudes via KLT relations. 

A more insightful "squaring" relationship between gauge theory and gravity amplitudes 
was conjectured recently by Bern, Carrasco and Johansson (BCJ) [2J. This relationship 
imphes a set of identities satisfied by gauge theory amplitudes, which were proven in [3]- 
[5]. BCJ conjectured that scattering amplitudes in gravity can be obtained from Yang- 
Mills scattering amplitudes, expressed in a suitable form, by a remarkably simple squaring 
procedure. In brief, the BCJ procedure involves expressing Yang-Mills amplitudes in terms of 
Feynman-like diagrams which only have cubic vertices. Associated with each cubic diagram 
is a colour factor and also a kinematic numerator. BCJ observed two remarkable things 
about these kinematic numerators. First, when a set of three colour factors satisfies a 
Jacobi identity, then the kinematic numerators can be chosen so that they also satisfy 
the same Jacobi identity. Second, the graviton scattering amphtude is obtained from the 

^ These three-point amplitudes are non-zero if the momenta of the particles are taken to be complex valued. 
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Yang-Mills amplitude by simply replacing the colour factors by the kinematic numerators; 
this "squaring" relationship has been proven in [H] assuming that appropriate gauge theory 
numerators exist. Explicit numerators have been described in Moreover, the relations 

are conjectured to hold off-shell ^U\. The BCJ relations have been discussed extensively 
in recent work, see [T0H2T] for example. The left-right factorization of the gravitational 
Feynman rules to all orders in perturbation theory has been recently demonstrated by Hohm 
[22]. 

It is especially intriguing that the kinematic numerators appearing in the BCJ relations 
satisfy Jacobi identities. Indeed, at four-point level this relation was found in the early 
1980s [SHIEI]. The fact that these Jacobi identities are satisfied at all points strongly suggests 
that there is a genuine infinite dimensional Lie algebra at work in the theory. Moreover, since 
the gravitational amplitudes are formed by replacing Yang-Mills colour factors with these 
kinematic numerators, this group seems to entirely determine the gravitational scattering 
amplitudes. It is worth emphasizing that the kinematic numerators which BCJ discuss are 
gauge dependent, so that the presence of this kinematic algebra seems to be manifest only 
in a class of gauges. In these gauges, the relationship between Yang-Mills theory and gravity 
is especially simple. 

In order to study this kinematic algebra in its simplest habitat, we consider the self- 
dual sectors of Yang- Mills theory and of gravity. The scattering amplitudes in these sectors 
vanish at tree level (with the exception of the three-point amplitude). However, a scattering 
amplitude is a rather perverse object to study in a classical field theory. It is more natural to 
study solutions of the field equations; indeed, tree-level scattering amplitudes can easily be 
obtained by taking a limit of a series expansion of an appropriate classical solution. We find 
that the kinematic algebra is manifest by a simple comparison of the equations governing self- 
dual Yang-Mills and gravity. The algebra is associated with area-preserving diffeomorphisms 
in a certain two-dimensional space. Our understanding of this algebra makes it completely 
trivial to see that the BCJ relations extend to the classical solutions of Yang-Mills theory 
and gravity which describe scattering of an arbitrary number of positive helicity particles, 
and measuring the resulting field at a point. There is an extensive literature on the hidden 
symmetries of self-dual gauge theory and gravity; see e.g. [^T] - Hi] . 

We further study the consequences of these results for the full Yang-Mills and gravita- 
tional theories. It is useful to consider these theories in the light-cone gauge, where they 
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can be formulated in terms of positive and negative helicity particles. Using a procedure 
inspired by Chalmers and Siegel [23], we are able to relate the MHV amplitudes to the self- 
dual sector, and to show that the corresponding BCJ relations are a result of the self-dual 
kinematic algebra. 

The structure of the paper is as follows. In section [IT] we open with a review of some 



background material. With this review in hand, we turn in sections III and IV to a discussion 
of the self-dual Yang-Mills and gravity theories, respectively. In these sections, we discuss 
the details of the kinematic algebra and BCJ relations for classical backgrounds. We discuss 
the implications of this algebra for the BCJ structure of MHV amplitudes of the full Yang- 



Mills theory in section |V| and of gravity in section |VI[ We close with some final remarks in 
section IVIII 



II. BACKGROUND MATERIAL 



Since this article deals with classical background fields applied to topics of interest in 
the recent literature on scattering amplitudes, it may be helpful to provide some review of 
our methods. We begin with a discussion of the BCJ relations, then provide an overview 
of the relationship between classical field solutions and tree amplitudes, and finally discuss 
the spinor-helicity method and a slight extension of this method which has been helpful in 
previous work on self-dual Yang-Mills theory and gravity. 



A. The BCJ relations 



Let us set the scene with a description of the BCJ relations [2]. The starting point is to 
express the M-point Yang-Mills amplitudes as a sum over the cubic diagrams, 

A. = /^-^Ef"' (1) 

where Cj are the colour factors, the denominators Dj are products of the Feynman propagator 
denominators appearing in the cubic graphs, and rij are kinematic numerators. The cubic 
diagrams can be thought of as diagrams showing how the Yang-Mills structure constants are 
sewn together to make the various appropriate colour factors for the scattering amplitudes; 
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at four points there are three diagrams (for the s, t and u channels), at five points there are 
15 diagrams and so on. Contact terms which appear in the Feynman diagram expansion of 
the amphtudes are assigned to appropriate cubic diagrams based on their colour structure; 
missing propagators in these contact terms are simply reintroduced by multiplying and 
dividing the contact term by such propagators. 

Given this decomposition, BCJ asserted that one can always choose numerators rij such 
that whenever the colour factors Cjj^,Cj^, and Cj^ of three diagrams ji,j2 and J3 are related 
by a Jacobi identity, then so are their kinematic numerators: 

Cji ± Cj2 ± = ^ Uj^ ± rij^ ± rij^ = 0, (2) 

where the appropriate signs depend on the definitions of the colour factors. Furthermore, the 
kinematic numerators can be chosen so that they enjoy the same antisymmetry properties 
as the colour factors under exchanging particle labels: 

Cj ^ -Cj rij -nj. (3) 

Taken together, these requirements indicate that the kinematic numerators nj have the same 
algebraic properties as the colour factors Cj. 

Given an expression for the Yang- Mills scattering amplitude in terms of numerators which 
satisfy the kinematic Jacobi identity, BCJ further noticed that the gravitational scattering 
amplitude is simply given by replacing the Yang-Mills colour factors by another copy of the 
kinematic numerator^, 

Mm = k''-'J2'^, (4) 

j 

where k is the gravitational coupling. 

It is worth remarking on the nature of the choices one has to make in defining appropriate 
numerators. Indeed, these numerators are not unique. For example, at four points the Yang- 

^ More general gravity theories can also be obtained using numerator factors from two different gauge 
theories [10] . 
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Mills scattering amplitude is written as 
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Since + Q + = 0, the amplitude can also be written in terms of new numerators 



where a is an entirely arbitrary function of the momenta and polarizations. Similar re- 
definitions are possible at higher points. BCJ described these redefinitions as generalized 
gauge transformations, since ordinary gauge transformations are a familiar way of moving 
terms between various Feynman diagrams. In the on-shell four-point case, the numerator 
redefinitions are such that the kinematic Jacobi identity, n'^ + + = 0, is maintained. 
However, as BCJ discuss, this is not true at higher points, even on shell. Thus, one must 
choose an appropriate gauge to find numerators which satisfy the kinematic Jacobi identity. 

B. Background fields as generating functions 

We know turn to reminding the reader of the connection between tree-level Feynman 
diagrams and classical field solutions [26]. In a few words, tree Feynman diagrams arise as 
the diagrammatic representation of terms in a perturbative series for the classical solution 
in the presence of a source, J . This can be understood from the functional approach in 
quantum field theory. Consider the generating functional 



where the fields are symbolically represented by 0. W[J] is the generating functional for 
connected correlation functions. Amplitudes are computed from the connected correlators 
by the LSZ procedure. Since our focus in this work is on tree level, we take the limit /i — > 0. 
In this limit, the path integral is dominated by solutions 0ci[^] of the equations of motion. 
Therefore, 






(7) 



w[j] = s[MJ]] + JMJ]- 



(8) 
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The n-point connected correlators are obtained by the nth functional differentiation of W[J] 
with respect to J. However, by the equations of motion, 

- +J)^^ + MJ] = MJ]- (9) 



5 J \ 6(j)c\ J 6 J 

Thus, the ra-point connected correlators are generated by n — 1 differentiations of the back- 
ground field created by the source. 

Much of the recent progress in our understanding of scattering amplitudes for massless 
particles has relied on the fact that the external legs of the amplitudes are on-shell, that 
is, their momenta p satisfy p^ = 0. There is a limit of these generating functionals for 
which all but one of the legs are on-shell. This technique has been used, for example, in 
the classic paper of Berends and Giele on recursion relations in gauge theory [27], as well 
as in a beautiful paper of Brown [2H] to sum tree graphs contributing to threshold particle 
production. Let us now illustrate how this technique works in the simple setting of scalar 
field theory. 

We will consider a massless cubic scalar field theory, with Lagrangian (including the 
source) 

C=^id<pf-^g<P' + J<p. (10) 
The classical equation of motion is simply 

d'(P+^gc^'-J = 0. (11) 

It is convenient to Fourier transform to momentum space; then the classical equation is 
written as 

k^<P{k) -\9 j ^Pi^P2 Kpi +P2-k) 0(pi)0(p2) = -Jik), (12) 
where we have introduced a short-hand notation 

m^i^Trmp). (13) 
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(1) 




FIG. 1: The first order correction (j)^^\k) consists of a single interaction between two particles 
sourced by J. 

We will solve this equation order-by-order in perturbation theory. Thus, we write 

oo 

0(fc) = ^0(")(fc), (14) 

n=0 

where (j)^'^\k) is of order g"^. At zeroth order, we must simply solve the free equation, so we 
learn that 

= (15) 
At next to leading order, counting powers of (7, we find that 0*^^^ satisfies the equation 

A;Vi)(A;) = ]^g j dpidp2Avi+V2 - fc) 0(°)(pi)0(°)(P2), (16) 



so that 



ik) = \l dp,dp, d{p, +P2-k)^^g (17) 

A J 1^ Pi P2 



Diagrammatically, this equation is shown in Fig. [TJ One can construct the three-point 
scattering amplitude from this expression by functionally differentiating twice with respect 
to the source J, and amputating the three lines. That is, 

■^3{pi,P2,P3) = -i lim lim lim pi pj pi —————— (j)'-^\-p3) = -igd{pi+p2+Pz)- (18) 

It is useful to note that the operation of differentiating with respect to sources and 
multiplying by is equivalent to differentiating with respect to the solution 0*^*^^ of the 
zeroth order equation. Indeed, one can then consider the limit 

J^k)^0, e^O, j(k)^-^j^O. (19) 
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In this limit, the equation of motion becomes 



5'0+^#' = O, (20) 

subject to the boundary condition that as (? — )■ 0, the complete solution tends to a non- 
trivial solution j{k) of the free equation (note that j{k) has support on fc^ = 0). Working in 
this limit has the advantage that the background field is built from the objects j{k) which 
have the interpretation of on-shell particles. The amplitude is given by 

A( ^ V 2 ^"-V("-^)(-Pn) .... 

A Pi, ...,Pn) = -^ hm —, r. (21) 

pi^O 5] [pi)5] {p2)...6j [pn-l) 

Note, however, that if we think of (j) not as a background field, but simply as a generating 



functional for tree amplitudes, then we do not need to take the limit (19). The difference is 
that, in the latter case, we can take all legs to be off-shell, while if (/> is a legitimate solution 
of the source-free field equation, then all but one legs (i.e. all legs representing sources) must 
be on-shell. 

Let us continue to one more order in perturbation theory. The next correction (j)^'^\k) 
satisfies 

k^<P^^\k)=g J (rp,ctp2d{p,+p2-k)<P^^\p,)<P^'\p2). (22) 
Inserting our expressions for (p^^^ and (f)^^\ we find 

If 11 

(p^'^\k) = - dpidp2dpz d{pi +P2+P3-k) 3{pi)j{p2)3{p3) 7^ 9 7 ^ ^ 9 (23) 

2y k^ [P2+P?,Y 

The Feynman diagram corresponding to this expression is shown in Fig. |2} The four-point 
amplitude is easily obtained by differentiation with respect to j and amputating the single 
remaining off-shell leg. 



C. Off-shell spinor-like variables 

In the next sections, we will apply these techniques to Yang-Mills theory and to gravity. 
In anticipation of these applications, it is useful to describe some slight extension of the 
spinor-helicity method. Let us begin in familiar territory. We define Pauli matrices 
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/ 1 \ 

/ (P2+P3)^ \ 

FIG. 2: The second order correction (j)^'^\k) consists of an interaction between two particles, 
creating a disturbance which propagates before scattering against a third particle. 

^-(;;)^ ^-(roO^ ^'-G-O' 

We also define = (cr°, — <?). These matrices satisfy the Clifford algebra 

a^d" + a^d^ = 27]^"". (25) 

Given a momentum p, one can consider the object paa = p-<^aa where a = 1,2 and a = 1,2. 
If = and p ^ 0, the matrix paa has rank 1. Thus we may write paa = AqAq, where A 
and A are two component spinors. Notice that A and A are not uniquely defined; indeed, 
new spinors \' = \/( and A' = (\ are just as good. We will exploit this freedom below. 

The only Lorentz invariant associated with a single momentum pi is pi = 0. However, 
given two momenta pi and p2, there are non-trivial invariants. We may form invariants from 
the spinors Ai and A2 associated with the two momenta. Lorentz transformations act on 
these spinors via a single copy of SU(2); therefore the object (12) = €"^XiaX2i3 is a Lorentz 
invariant. Similarly, the object [12] = e^^Ai^Ag^ is an invariant. All other Lorentz invariants 
are built from these. In particular {pi +^2)^ = S12 = (12)[12]. 

It will be convenient for us to define analogues of the spinorial products for off-shell 
momenta. In doing so, we follow the notation of |29l |30] . We exploit the freedom to rescale 
the spinor A associated with an on-shell momentum p to write A = {Q, 1). Note that A 
is now dimensionless; correspondingly, A has (mass) dimension 1. We then compute that 
>^ = {Pw,Pu) where 

u = t — z, V = t + z, w = X + iy. (26) 
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The spinor products are then given by 



[12] -> X{pi,p2) = PlwP2u - PluP2w, (27) 

(12) ^ Qipi,P2) = Q{pi) - Q{P2). (28) 

We will use the notation X{pi,p2) and Q{pi,P2) to emphasize that the spinor product is 
taken with this unusual rescaling of the spinors. The benefit is that the definition of X 
extends to arbitrary off-shell momenta. We shall have more to say about the freedom to 
rescale spinors in the later sections of this paper. 



III. SELF-DUAL YANG-MILLS THEORY 



In this section, we begin our study of classical background fields. Since our goal is to 
investigate any possible BCJ-like structure of the fields, we choose to study the simplest non- 
trivial fields available. These are the self-dual solutions. Firstly, we consider the self-dual 
Yang-Mills (SDYM) equations in Minkowski spacetime, 

Pfiu /ii/paP^ • (29) 

The gauge field is necessarily complexified, and the physical interpretation is that it is a 
configuration of positive helicity waves. Our set-up follows Bardeen and Cangemi 



see also [M] . We work with the coordinates ( 26 ) , such that the metric is given by 



ds'^ = dudv — dw dw. (30) 



We choose the light-cone gauge, where Au = 0. The self-dual equations (29) then imply 



A^ = 0, A, = -^d^<!>, A^ = -^du^, (31) 



^ In our conventions, the field strengh is F^^^ — S^A^ — d^A^ — ig[A^, A^], and the structure constants of 
the Yang-MiUs Lie algebra are defined by [r'^,T^] = if^'^T". 
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and also an equation of motion for the Lie algebra- valued scalar field $, 

+ = 0, (32) 

where (9^ = 4 — dwdw) is the wave operator. Thus, our problem reduces to studying a 
scalar equation with a cubic coupling. 



Following the discussion in the last section, let us now solve the scalar equation (32) with 
the boundary condition that, when (7 — )■ 0, $(a;) — ?■ j{x). In momentum space, we can write 
this equation as 

^\k) = \gj ^p,^P2 ^'''\r'''^ <^''^ipi)^''iP2), (33) 
where we have defined 

Fp.p,'' = d{pi+p2-k)X{pi,p2). (34) 
We shall use an integral Einstein convention for the contraction of the indices of Fp^pJ^, 



F " F ^ 



j dq d{pi + q - k) X{pi,q) d{p2 + Ps - q) X{p2,Pz) 
= 3(pi +P2+P3- k) X{pi,p2+P3)X{p2,P3)- (35) 

Moreover, we can lower and raise indices using 

= d{p + q) = 6pg, such that dpgd'^'' = 5p^ = d{p - k). (36) 
It is straightforward to see that F^^p^ps = Fp^p^p^ is totally antisymmetric, e.g. 

PP.P2P. ^ +p3)X(p,, -pi -P3) = +p3)X(pi,p3) = -F™^ (37) 

Our notation is designed to emphasize the fact that the coefficients F^iP2P3 have the same 
algebraic properties as the structure constants f°'^'^. However, we will learn below that 
the significance is deeper, and that F^iP2P3 qjq^ fact, structure constants for a certain 
infinite-dimensional Lie algebra. 



Equipped with this formalism, we solve the equation of motion for (33) iteratively, as 
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J 




$(3) 



J 



FIG. 3: The third order correction ^^^\k), in (41), consists of two different terms. The first one 
corresponds to the Feynman diagram on the left, and the second one to the diagram on the right. 



a series expansion in the coupling constant g. The first few terms are 



^^'^%k) = -g j dp,dp2^^ j'HpOj'Hps 



rp k rp q fbica £h2ho^c 



63 1 



As one would expect, different terms appear at higher orders. 



(38) 
(39) 

(40) 



$(3)«(fc) 



F 91 F 92 f feicia f fe2C2Ci fhzb4,C2 

J- nnnn J- noriA J J J 



P3P4 



k'^{p2 +P3 +P4)^(P3 +. 



+ 



+ 



F F 91 F 92 fciC2a fbib2Ci f 6364C2 
9192 PiP2 P3P4 J I J 

Ak'^{pi +P2Y{P'3+PiY 



f'{Pl)f'{P2)f'{p3)j''{p,). (41) 



The two terms correspond to two types of Feynman diagrams, represented in Fig. [3] 

The way in which the kinematic factors F^^^^^^ mirror the commutators is a consequence 



of equation (33). Indeed, the colour/kinematics correspondence is simply 



abc 



(42) 



In the same way that each vertex is dressed by a /"'"^ factor, it is also dressed by a kinematic 
ppiP2P3 factor. Finally, the kinematic Jacobi identity is given by 



F IF ^ + F F ^ + F F — n 

J- piP2 ^ psq ' ^ P2P3 ^ Pl9 ' ^ P3Pl ^ P29 ~ ^1 



(43) 
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which is a consequence 



X{h, k2)X{h, h + k2) + X{k2, ks)X{k,, k2 + ks) + X{ks, k,)X{k2, ks + k,) = 0. (44) 

The relation that we found here is exactly of the type that BCJ pointed out for Yang- 
Mills scattering amplitudes [2J, but now in the context of classical background solutions. 
The obvious question now is: what is the algebra whose structure constants are the F^^^^pa 
factors? The answer appears naturally once we look at the gravity case. 

For completeness, let us remark that there is a simple expression valid to all orders [20], 

= (ig)^ j dpid:p2 . . . ipn+l + ^2 + • • • + Pn+l - k) 

X j{Pl)j{P2) ■ ■ -jiPn+l) Q{Pl,P2y^Q{p2,P3y^ ■ ■ ■ Q{Pn,Pn+lY^- (45) 

However, this form is not convenient for our purposes, since the Feynman diagram expansion 
is not explicit. Moreover, unlike the preceding discussion, this expression is only valid if j{k) 
has support on k'^ = 0, i.e. if the legs representing sources correspond to on-shell particles. 



IV. SELF-DUAL GRAVITY 



We consider now self-dual gravity (SDG) with Lorentzian signature, using an approach 
analogous to the one we used for SDYM. Previously, Mason and Skinner ^2] have computed 
the MHV amplitudes in gravity by perturbing the classical, self-dual background field, in 
work analogous to that of Bardeen and Cangemi [29l |30] on self-dual Yang- Mills theory. Our 
focus here will be on the background field itself rather than perturbations around it. The 
self-dual equations are 



Inspired by the gauge field (31), we try as a solution the metric 

g^lu = Vfj-u + t^h^^, (47) 



On-shell, this identity can be understood as a special case of the Schouten identity. 
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where k is the gravitational couphng, ?7^j, is the Minkowski metric (30 ), and the non- vanishing 
components of h^y are 



hr. 



4 



1 



(48) 



The SDG equations (46) then imply that the scalar field (p obeys 



d'<P + K{{dl<f)) (9» - {d^du<Pf) = 0, 



(49) 



where denotes the Minkowski space wave operator. It turns out that this equation was 
first obtained by Plebanski [33], and that it allows for the most general solution of SDG. 



as 



The resemblance between SDYM and SDG becomes even more striking if we rewrite (49) 

d^(j) + K{d^<j),du(l)} = 0, (50) 



which should be compared to (32). We introduced here the Poisson bracket 



{f,9} = {dwf) {dug) - {duf) {du,g), 



(51) 



from which we construct the Poisson algebra 



(52) 



This is the kinematic algebra that we were looking for in the last section. It is the Poisson 
version of the algebra of area-preserving diffeomorphisms of w and u. To see this, consider 
a diffeomorphism w — w'{w,u), u — ?■ u'{w,u). This transformation preserves the Poisson 



bracket (51) if and only if it has a unit Jacobian, i.e. it is area-preserving. The infinitesimal 



generators of the diffeomorphisms are 



Lk — 6 ( kyjdu + k^d^})- 



(53) 



The possible contributions to the right-hand-side of ( 49 1 can be absorbed by a redefinition of 
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The Lie algebra is 

[Lp^,Lp^] = iX{pi,p2)Lp^+p, = iFp^pJ'Lk, (54) 



and the Jacobi identity was given in (43). 

We have seen that this algebra is the kinematic analogue of the Yang-Mills Lie algebra. 
It is interesting to note that there is a correspondence between the Lie algebra of area- 
preserving diffeomorphisms of S"^ and the Lie algebra of the generators of SU(A^) in the 
planar limit — )■ oo, in the sense that there exists an appropriate basis such that the 
structure constants are the same 



Let us now proceed, as in the SDYM case, to solve this gravitational equation (50). In 
momentum space, the equation can be written as 

m = ^kJ dprdp, ^^P^^P^^^p^P'^' 0(pi)0(p2), (55) 

We again take (p^^^k) = j{k) to have support on the light cone. To the first few orders in 
the solution is 



0W(fc)=j(A:), (56) 
1 [ ^ ^ X{pi,p2) 

-kJ dpidp2 p 

j dp.dp^dp,'^ P(p.+P3)^ 



«(A:) = 1^1 dp,dp2 ^^^^^^^i^j{pMP2), (57) 
0(^)(A:) = ln' [ dp,dp2dp, ^^^''^^fr;'^^^^:r^^""% -(Pi)j(P2)j(P3). (58) 



Let us explore the relationship between the expressions (33 40) and (55 58). Indeed, 
it is clear that one can deduce the gravitational expressions from the Yang-Mills cases by 
replacing the SU(N) structure constants /"''^ by appropriate factors of X. These factors of 
X are, in turn, related to the structure constants F of the kinematic algebra. However, the 
relationship is not given by / — > F because this would involve squaring a delta function. One 
algorithm for deducing the gravitational expressions from the Yang-Mills formulae involves 
extracting the overall momentum conserving delta function, and then following the BCJ 
procedure of identifying a kinematic numerator which is to be squared. Let us illustrate 
this at the level of the second corrections, and Beginning with the Yang-Mills 
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formula, we have 



$(2)a(^^^| =]^g^ j dpidp2dp^ d{pi +P2+P3-k) 

fbica fb2hsc 

f'iPl)f'iP2)f'iP3)- (59) 



The corresponding off-shell tree-level correlator - we will refer to this object as an off-shell 
"amplitude" - is obtained by differentiating with respect to the and amputating the 

final leg. The expression now becomes 

A'^' {PuP2,P3,-k) = - g — ^ — ^{pi+P2+P3-k)+- ■ ■ , 60 

2 {P2+P3>^ 

where the dots indicate two other channels. The BCJ procedure now identifies the kinematic 
numerator as 

nt = X{pi,p2+P3)X{p2,P3) (61) 

This is the object which must replace the colour factor fbicajb2bzc deduce the gravita- 



tional perturbation, as one easily checks. Moreover, as expression (44) shows, these objects 
satisfy a Jacobi identity. Of course, we could have deduced this numerator directly from 
the background field perturbation; we have introduced the off-shell "amplitude" merely to 
illustrate that the procedure is equivalent to the BCJ procedure. 

V. YANG-MILLS THEORY 

We have seen that there is a remarkably close connection between the equations governing 
self-dual background fields in Yang-Mills theory and in gravity. We will now consider the 
full Yang-Mills theory in the light of the kinematic algebra present in the self-dual sector. 

There is one non-trivial scattering amplitude in the self-dual theory, which is the three- 
point amplitude (evaluated for complex momenta). In fact, this three-point amplitude is the 
same as the amplitude in the full theory. Let us briefly explain why this is so. In Yang-Mills 
theory, the three-point + H — scattering amplitude is given by 

(pi,P2,P3) = -^g r''d{p,+p2+P3). (62) 
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where f"'^'^ = It is convenient to rescale the spinors as described in section |ll] to 

make contact with the invariants X{pi,pj). In terms of a scaling factor for z = 1, 2, 3, we 
have 

[12] = ClC2X(pi,p2), [23] = C2C3X{P2,P3), [31] = C3ClX{p3,Pl), (63) 

so that the three-point amplitude can be written as 

At{pi,P2, P3) = -^9 ^ ^^^\'^f J "'' API +P2+ P3) (64) 

Cs X{p2,P3)X{p3,pi) 

= -^gS^X{p,,p2)r'"'^{p,+p2+P3) (65) 
^-tgF,,p,p,r'"^, (66) 

where in the last line we redefined our basis of polarization vectors to absorb the factors Q. 
Then this Yang-Mills amplitude is given by the same expression as in the self-dual sector. 
It is important to note that, while we started from a spinor-helicity on-shell expression, the 



final result (66) can be obtained off-shell. 

Similar comments hold for the h Yang-Mills amplitude; a basis of polarization vectors 

can be found so that the YM amplitude is given by the three-point amplitude in the anti- 
self-dual sector. It is convenient to define 

-^pip/ = 5(pi +P2- k) X{pi,p2), where X{pi,p2) = piwP2u - P\uP2w- (67) 

Then the anti-self-dual three-point amplitude (and the V YM amplitude) is simply 

(pi,P2,P3) = 2^7X(pi,p2)r''^5(pi +P2 +P3) = (es) 

where -Fpipapg = X(j)\,p2)t{j)\ + p2 + Ps) are the structure constants for the anti-self-dual 
kinematic algebra. Note that, if p\ and p2 are on-shell, then the quantities X(pi,p2) and 
X(pi,p2) are related by 

X{px,P'2)X{px,P2) = -^PluP2uSl2- (69) 

In the gravitational case, it is also true that a polarization basis can be chosen so that 
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the three-point functions in the full theory are identical to the three-point functions we have 
discussed in the self-dual sector. 



A. Chalmers-Siegel Theory and the Four-Point Amplitude 

We have seen that the kinematic algebra is present in the full theory, at least at the level 
of the three-point amplitudes. However, it seems that two algebras appear: one for the 



MHV three-point amplitude, and another for the MHV case. To study how these algebras 
interact, it is convenient to turn to the Chalmers-Siegel formulation of Yang-Mills theory in 
light-cone gauge [25]. Their action describes Yang-Mills theory as a theory of interacting 
positive and negative helicity degrees of freedom, and therefore is very helpful to connect 



the full theory with the discussion of Section III We will see here how the colour/kinematics 



duality for the self-dual (or anti-self-dual) sector underlies the BCJ relations in the MHV 
sector of Yang-Mills theory. 

As before, we take the light-cone gauge 

A^ = (0,v4„AA), (70) 



where we again work with the coordinates (26), and denote A = A^,, A = A^,, to avoid the 



proliferation of indices. Positive helicity particles correspond to A = 0, and negative helicity 
particles to A = 0. The great simplification of the light-cone gauge is that now A^ appears 
quadratically in the Yang-Mills Lagrangian and can be integrated out. We are left with 

£ = tr{ }■ (71) 

The equations of motion are 

ld'A = - ig [d^A, A] + ig d^A] - ig [^A + ^A, duA] 

+ g'[A, hA,d.A]]+g'[d.A, hlAduA] + [A,d.A])] (72) 



and the conjugate equation for A. One can easily check that equation (32) is recovered for 
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the choice A = 0, A = — 9„$/4. In momentum space, we have 

h W ^ fbib2a 
"'u ^ V1P2 J 



J PluP2u 1^ 

Plu PpiP2 . 

P2uku k"^ 



+ 4 ^ '-^^ A'^ {pM'^ iP2) I 



+ Ag' [ dp,dp,dp,d{p,+p,+p, - k) +Pl"P3n f ^.,^^^)^.,^^^)^63(p^)_ 

J {P2u+P3u)^ 

(73) 

There are now two more vertices with respect to the self-dual case: the cubic vertex from 
the anti-self-dual theory and the four-point vertex. 

We proceed as in the past sections to construct the solution perturbatively, starting from 

A(o)<'{k) = 7]^{k), A^^'^'^ik) = f^^ik). (74) 

We can then obtain the tree-level -|- H (off-shell) amplitude 

.4(1*2-3-4-) = (75) 

The on-shell amplitude is computed by taking the limit pj — )■ 0, for i = 1,2, 3, 4. There are 
several Feynman diagrams contributing to this amplitude, including the four-point vertex 
diagrams. We choose a gauge by taking p4 to be a reference leg in the sense of Chalmers 
and Siegel [25]. Thus we take the limit p^u 0,^ so that the four-point vertex is eliminated. 

^(1+2+3-4-) = g' sUpi)eUP2)e^iP3KM — x 

Piu PluP2u 

7 pi fbib2C fbsbiC rp q p fb2bzc fb-ib^c rp g p f bs^ic f 62^4 

PlP2g ^ P3P4 J J 1^ ^ P2PS ^ PlPiQ J J |_ ^ PSPl ^ P2P4,q J J 

S U t 

where we introduced the Mandelstam variables, 



(76) 



s = (Pl+P2)^ ^ = (Pl+P3)^ u = {pi+pif. (77) 



We consider complex momenta, and thus this condition does not imply that p4 is on-shell. 
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There is no divergence in the amphtude as p^u — ^ since it is cancelled by eyj{p4) oc p4^u- 



The form of the amplitude (76) indicates that the BCJ kinematic numerators are 



Us = aX{pi,p2)X{p3,pi), nt = aX{p2,P3)X{pi,p^), = aX{p3,pi) X{p2,Pi), (78) 

where a is a proportionality factor. Because we now have both X and X, it may not seem 
that 

+ nt + riu = 0, (79) 

but recall that we set p4u = 0, so that 



Us = a' X{pi,p2) X{p3,p4), Tit = a' X{p2,P3) X{pi,p4), Uu = a' X{p3,pi) X{p2,P4), 



(80) 



where a' = ap^ui/p^w- Then the relation (79) is nothing but the Jacobi identity (44) of the 



self-dual theory. We stress that all external legs can be taken to be off-shell. 
In order to go on-shell, we can take p^u = Paw = and write 



ns = -ap4iijX{pi,p2)p3u, nt = -ap4u,X{p2,P3)piu, = -ap^a X{p3,pi) p2u, (81) 



so that ( 79 ) now becomes 



X{pi,P2)p3u + X{p2,P3)Plu + X{p3,pi)p2u = 0. 



(82) 



Ref. [25] details the procedure by which the amplitude can be put in the standard spinor- 
helicity form, so we will not address this. 



B. Higher-point MHV amplitudes 

In our study of the four-point Yang-Mills amplitude, we found it helpful to make the 
reference leg choice p^u 0. Any reference leg must be attached to a three-point vertex 
of the same self-duality; that is, a negative helicity reference leg must be attached to a 

H vertex. In the case of the four-point function, this choice of reference leads to the 

off-shell amplitude being written in terms of products of a single H vertex and a single 

+ H — vertex. Moreover, with this choice both vertices are proportional to the structure 
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constants of the self-dual kinematic algebra. We shall again exploit this simplification to 
clarify the relationship between the Jacobi identities satisfied by the MHV amplitudes and 
our kinematic algebra. 

We will show that all the Feynman diagrams contributing to an MHV n-point amplitude, 
containing n — 2 particles of positive helicity and 2 particles of negative helicity, contain 
only three-point vertices. By choosing a negative helicity line to be the reference leg, we 

ensure that the vertex attached to that leg is a h vertex, while all the other vertices 

are + -\ — vertices. To see this, let n_|_ be the number of + H — three-point vertices in a 

given Feynman diagram, while n_ is the number of 1- vertices. Similarly, let 71,4 be 

the number of four-point vertices. We suppose there are / internal lines at n points. It is 
straightforward to establish a number of identities relating these quantities. Counting the 
powers of the coupling g, we find that 

n-2 = n+ + n^ + 2n4. (83) 

Meanwhile, the positive helicity ends of internal or external lines must terminate in the 
vertices, and similarly for the negative helicity lines: 

n-2 + I = 2n+ + n_ + 2n4, (84) 
2 + / = n+ + 2r2_ + 2^4. (85) 

Manipulation of these identities leads to the conclusion that there is a strong constraint on 
the sum of the number of h vertices and four-point vertices: 

n_+n4 = l. (86) 

Since we choose one negative helicity line as a reference, it follows that n_ = 1 and so 
?T,4 = 0. Therefore, the MHV amplitude is computed from cubic diagrams. As an example, 
we present in Fig. |4] a diagram contributing to the 6-point MHV amplitude. The negative 
helicity line represented in bold is the reference leg, which is connected to the only h 
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4+ 



FIG. 4: Example of a diagram contributing to the 6-point MHV amplitude. The thick line repre- 
sents the reference leg; the vertex attached to it is of the type h, while the other three vertices 

are of the type + H — . 



vertex. The contribution from this diagram is proportional to 

5C3 fCiC2C3 

(87) 



TP qi TP 92 p 93 P fbib^ci fb2bzC2 fbib^cs fciC2C3 

^ PlP6 ^ P2P3 ^ P4P5 ^ 919293 J J J J 



{pi + PfsY{.P2 + pzY{,Pi + p^y 

Just as in the four-point case, the choice of reference hne attached to the single h 

vertex in the amplitude means that the quantity X appearing in the vertex is proportional 
to X. Expressed in this way, the Jacobi relations satisfied by the numerators are nothing 
but the Jacobi relations of the self-dual kinematic algebra. 

VI. GRAVITY 

In this section, we will see how the kinematic algebra arises for gravity, and confirm that 
it indeed corresponds to the BCJ "square" of the Yang-Mills case. We need a formulation of 
the theory in terms of positive and negative helicity particles, analogous to the Chalmers- 



Siegel action, in order to make a connection with Section IV Fortunately, such a Lagrangian 
was obtained by Ananth et al |35]; see also [36] . 

The Lagrangian is obtained from the Einstein- Hilbert Lagrangian in the following way 
(see Appendix C of [35j for a detailed derivation). Take the light-cone gauge for the metric, 

Quu = gui = 0, (88) 
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where i = 1, 2 denotes the x and y Cartesian coordinates. Parameterize the remaining 
components as 

9u. = \e^'\ ^.• = -e^/S„ (89) 

where is a symmetric 2x2 matrix of unit determinant. Then the Einstein constraints 
Ruu = Rui = allow us to write the Lagrangian in terms of 7jj only. The perturbative expan- 
sion in terms of positive helicity gravitons h and negative helicity gravitons h is performed 
by taking 

, „, I I h + h i(h-h)\ 

^ = (e-^).. , H=—\ ^ ' \ . (90 

' ^ V2 y^{h-h) -{h + h)) ^ ^ 

Expanding in the coupling k, we obtain the Lagrangian^ 

C=-hd^h + Khdf, -h$^h+(^h]\+Khdf, -h%h+(^h ' 



1 - d d- - 1 

— {duhduh) -^{duhdji) + — 



{duh duh) [duh duh) + — [duh duh) {dyjdyjh duh + duh dwd^h) 



- 7^ {duh duh) 2 du,du,h h + 2h d^d^h + 9 d^h d^h + du,h d^h ^^^hduh - duh h 

' \ Ou Ou J 

- 2^ (2 dyjh duh + h dwdji — d^duh h) h du,h — 2^ (2 duh d^h + du,duh h — h du,duh) dyjh h 

Ou Ou 

- ^ (2 du,h duh + h du,duh - du^duh h) ^ (2 duh dji + du,duh h-h du,duh) 

Ou Ou 



h h ( dujdujh h + h d^d^h + 2 du,h d^h + 3 ^ ^ h duh + 3 duh ^ ^ h 

du du 



+ 0{k^). (91) 



The self-dual equation (50) is easily recovered from the general equations of motion by 
setting h = 0, h = -dl(l)/A. 



^ As described in |35j . a field redefinition which removes occurences of 9„ from the four-point interaction 
has been performed; this is irrelevant for our purposes. 
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In momentum space, the equations of motion take the form 

J ^ PluP2u 



k 



kl X{p,,p,)F, 

PluP2u 

, , Pin ^iPuP2)Fp,p^ - 

+ 4 p Kp,)h{p2] 

+ four-point vertex + 0{k^). (92) 



We aheady see the squaring of the kinematic factors in the cubic vertices, when comparing 



with the analogous expression in Yang- Mills theory, equation (73). More explicitly, we can 

determine the tree-level + H (off-shell) amplitude. We saw in the Yang-Mills case (76) 

that the reference leg choice — t- was very useful. The same procedure leads now to 

16 

A^(l+2+3"4^) = e+^{pi)e+^{p2)e^i;,{p3)€^^{pi) — d{pi + p2 + Ps + Pa) 

Pau PiuP2u 

' Xip,,p2yxips,p,f ^ X{p,,p2fX{p,,p,f ^ Xip,,p2fX{ps,p,f 
s u t 



The X vertex associated with the reference leg is proportional to X, as we saw before. The 
amplitude is finite, due to the contribution from the polarization vectors, which are taken 
to be the squares of the Yang-Mills polarization vectors. 



e 



{k)=e^{k)e^{k). (94) 



Notice that it is crucial that the four-point vertex vanishes in the limit — ?■ 0. This can 



be seen directly from the Lagrangian (91): never acts on a single particle {h or h) in 

the contribution. 

It is no surprise that the BCJ squaring ^ holds for this MHV amplitude. In fact, under 
a certain assumption, we can show that it holds for any MHV amplitude. The argument is 



the same as the one used in Section VB to show that the four-point vertex is irrelevant for 
MHV amplitudes, if we use the reference leg gauge. Again, we consider an n-point MHV 

amplitude + ■ ■ ■ H , and denote the number of internal lines by /, the number of + H — 

vertices by and the number of H vertices by n_. The additional vertices are labeled 

by the number of incoming positive helicity lines o"+ and the number of incoming negative 
helicity lines cr_; the number of such vertices is denoted by n„_^„_. These additional vertices 
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all satisfy 0"+, cr_ > 2, which can be explicitly checked by extending the Lagrangian above 



to higher orders in k. Then, the expressions analogous to (83 85) are 

n - 2 = n+ + n_ + ^ (a+ + a_ - 2) n„_^_a_ , (95) 

cr-j-,0"_ 

n - 2 + / = 2n+ + ?7._ + ^ (T+ ?7.^^^_ , (96) 



2 + J = n+ + 2n_ + ^ (T_ (97) 



a+,a- 



from which we conclude that 



+ 5^(a_-l)n.^._ = l. (98) 



(J -I- ,(T 



We assume that, as in Yang-Mills theory, the reference leg must be attached to a cubic 
vertex. We choose a negative helicity particle as our reference leg, so that n_ = 1. It follows 
that the MHV amplitudes can be computed only from cubic vertices. For our assumption 
to be true, it is required that no higher-order term in the Lagrangian contains 1/9^ acting 
on a single particle. Ref. [37] presents the contribution to the Lagrangian, and this is 
indeed the case. We conjecture that this is valid to all orders. 



VII. CONCLUSIONS 



In the early sections of the paper, we described how a study of the self-dual sectors of 
Yang-Mills theory and of gravity sheds light on the BCJ relations. Cubic Feynman diagrams 
appear quite naturally in the self-dual sectors, and it is easy to see that the numerators of 
these diagrams have the properties which BCJ described. This occurs for a very simple 
reason: the numerators are built out of structure constants of an infinite dimensional Lie 
algebra. In the case of gravity, the presence of an infinite dimensional algebra seems entirely 
reasonable, and it is pleasing to see such an algebra play the role in gravity that the finite 
dimensional colour algebra plays in Yang-Mills theory. Indeed, Mason and Newman [38j (see 
also [3^ ) previously proposed a very similar relationship between gauge theory and gravity. 
But it is quite remarkable that the KLT relations somehow work in reverse, so that this 
algebra is also present in the Yang-Mills theory. One of the main results of this paper is that 
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the Yang-Mills cubic vertices are dressed by the straight product of the structure constants 
of the colour algebra and the kinematic algebra. 

To make the kinematic algebra manifest in the case of MHV amplitudes of Yang-Mills 
theory, it was useful for us to work in the Chalmers-Siegel formulation of the theory. By 
choosing a special gauge, we removed the four-point vertex, so that the MHV amplitudes are 
computed off-shell by cubic diagrams. This is reminiscent of the BCFW computation of the 
MHV amplitudes. The close connection between BCFW and the Chalmers-Siegel action has 
already been pointed out by Vaman and Yao [17]. It seems that to clarify the structure of 
the kinematic algebra beyond MHV amplitudes will require understanding how to formulate 
these amplitudes in terms of cubic diagrams. Such a formulation must be related to the pure- 
spinor calculation of the BCJ numerators described by Mafra, Schlotterer and Stieberger 
[9] . Similarly, it would be interesting to see the supersymmetric version of this algebra. The 
work of Broedel and Kallosh |10] on the light-cone formulation of = 4 SYM and A/" = 8 
supergravity would be a good starting point. 

On the gravitational side, we have seen how to compute MHV amplitudes from cubic 
diagrams in light-cone gauge, up to a mild assumption about the behaviour of the higher 
terms in the gravitational light-cone Lagrangian. The fact that the BCJ relations hold (at 
least on-shell) for all graviton amplitudes suggests that there is some cubic theory which can 
be used to compute these amplitudes. Finding this theory would be a great step forward in 
our understanding of both gauge theory and gravity. 
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